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plications derived from recent microwave surveys, with emphasis to those coming from 
the Planck mission. We critically discuss the impact of systematics effects and the role 
of methods to separate the cosmic microwave background signal from the astrophysical 
emissions and each different astrophysical component from the others. We then review 
of the state of the art in diffuse emissions, extragalactic sources, cosmic infrared back- 
ground, and galaxy clusters, addressing the information they provide to our global view of 
the cosmic structure evolution and for some crucial physical parameters, as the neutrino 
mass. Finally, we present three different kinds of scientific perspectives for fundamental 
physics and cosmology offered by the analysis of on-going and future cosmic microwave 
background projects at different angular scales dedicated to anisotropics in total intensity 
and polarization and to absolute temperature. 

Keywords: cosmology, cosmic background radiation, galaxy clusters, active galaxies, pri- 
mordial galaxies, Milky Way, Zodiacal Light. 

PACS numbers: 98.80.-k, 98.70.Vc, 98.65.-r, 98.54.-h, 98.54.Kt, 98.35.-a, 96.50.Dj. 
1. Introduction 

Since its discovery, the cosmic microwave background (CMB) represents a crucial 
probe for our general view of the Universe and the understanding of key aspects in 
cosmology and fundamental physics. Furthermore, microwave surveys are becoming 
more and more relevant for the comprehension of the physical and evolutionary 
properties of astrophysical structures at different cosmic epochs, from galactic to 
cosmological scales. Following the very important results from balloon-borne exper- 
iments, the NASA COBE and WMAP satellites, and recent ground-based projects, 
covering together a wide multipole range, the available and forthcoming data prod- 
ucts from the Planck missiorQwill have a strong impact in these fields in the coming 
decade^] Planck instruments are in fact the most sensitive microwave receivers ever 
launched in space. Their sensitivity calls for a comparable level of systematic effect 
control, one of the main drivers in satellite and instrument design and currently key 
in data reduction and interpretation^^, a topic addressed in Section [2| Similarly, 
for high sensitive microwave observations, the accuracy in the recovery of the CMB 
properties largely relies on the capability to disentangle the cosmological signal from 
the astrophysical emissions (the so-called foregrounds) , as discussed in Section [3] 
Waiting for the Planck cosmological results, we focus here on currently available 
astrophysical discoveries based on the first Planck surveys possibly complemented 
by other sets of observations carried out at similar wavelengths and combined with 
surveys in other frequency domains, in a multifrequency approach. Section [4] is de- 
voted to the diffuse emissions coming from the Solar System and the Galaxy, on the 



a Planck is a project of the European Space Agency - ESA - with instruments provided by two 
scientific Consortia funded by ESA member states (in particular the lead countries: France and 
Italy) with contributions from NASA (USA), and telescope reflectors provided in a collaboration 
between ESA and a scientific Consortium led and funded by Denmark. 

b This paper is based largely on the Planck Early Release Compact Source Catalogue and pub- 
licly available publications by ESA and the Planck Collaboration, for what concerns the related 
aspects. Any material presented here that is not already described in Planck Collaboration papers 
represents the views of the authors and not necessarily those of the Planck Collaboration. 
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physical processes operating at various Galactic scales, and on their manifestation 
at microwave wavelengths. In Section [5] we discuss the main properties derived on 
extragalactic sources at different cosmic distances, including the main evidencies 
regarding the cosmic infrared background. Section [6] is dedicated to the recent ob- 
servational results on galaxy clusters and on their scaling relations, with the wealth 
of information they provide on baryon physics and their use as a probe for neu- 
trino mass estimation. Finally, in Section[7]we review three different topics in CMB 
cosmology, possibly linked to fundamental physics, that will be addressed respec- 
tively by the forthcoming results from the Planck mission, by future high resolution 
ground-based experiments, and by the next generation of CMB spectrum projects. 

2. Control, assessment and removal of systematic effects in Planck 

Planck orbits around the L2 Lagrangian point and scans the sky spinning at 1 rpm in 
almost great circles with its Gregorian dual-reflector telescope pointing at 85° from 
the spin axidSE. In the telescope focal plane the microwave photons are collected by 
two wide-band receiver arrays spanning a frequency interval ranging from ~30 GHz 
to ^857 GHz. The Low Frequency Instrument (LFI), is a coherent differential array 
based on 20 K InP HEMTQ amplifiers currently working in three bands centered at 
approximately 30, 44 and 70GHz^. The High Frequency Instrument is an array of 
bolometers cooled to 0.1K operating at six frequency bands centered at 100, 143, 
217, 353, 545 and 857 GHz 8 . Planck (full width half maximum, FWHM) resolution 
ranges from 33.3' to 4.3' going from 30 GHz to 857 GHz, and its final sensitivity 
per (FWHM 2 resolution element is in the range of ~ 2—14 [iK/K in terms of 5T/T 
for frequencies v < 353GHz. The life of Planck largely exceeded the early plan. Five 
all-sky surveys has been accumulated with HFI, while LFI is planned to operate up 
to about the end of Summer 2013, so completing eight all-sky surveys. 

In Table [l] we list the main systematic effects in Planck according to their source 
and provide few notes about their control and residual impact on science. 

In the LFI we have generated timelines, maps and power spectra of thermal 
effects and 1-Hz spikes using in-flight scientific and housekeeping data coupled with 
transfer functions measured during ground tests. Table [2] reports the peak-to-peak 
and rms effect on full-sky temperature maps of these effects, while in Fig. [I] we show 
their expected temperature angular power spectra after component separation. The 
effect of component separation (see next section) has been reproduced by mixing the 
systematic effects maps using the same mixing matrix used to extract the Planck 
CMB map (not reported in this paper) from the individual frequency maps. 

Our current analysis confirms that the level of systematic effects rejection is 
in line with pre-launch expectations and will allow full exploitation of the science 
encoded in the CMB signal. 



c Indium Phosphide High Electron Mobility Transistor 
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Table 1. Main systematic effects in Planck 
Category Effect Notes 

OpticMm 

Side-lobes Galaxy and CMB dipole pickup by main and sub- 

reflector spillovers. Negligible effect on tempera- 
ture maps, needs to be removed at low frequency 
for polarization analysis. 

Detectors 

Cosmic ray hits Affect bolometric detectors. Removed from time- 

lines via template fitting. 

1/f noise Affects radiometric and bolometric detectors. In 

the LFI the 1/f contribution is limited to max 
3% by differential measurement strategy and de- 
striping algorithmic. 

Bandpass mismatches' Affects primarily radiometric detectors. Negligible 
impact on temperature. Corrected in polarisation 
at map level exploiting polarized source measure- 
ments (Crab Nebula). 

Electronics 

1-Hz spikes Affects LFI data. Removed from timelines by tem- 

plate fitting. 

Thermal^ 

300 K fluctuations .... In principle affect both instruments. Inherent 
hardware stability is compliant with scientific re- 
quirements. 

20 K fluctuations Affect mainly LFI. Inherent hardware stability is 

compliant with scientific requirements. 
4 K fluctuations Affect both instruments. Inherent hardware stabil- 

ity is compliant with scientific requirements, 
tin the LFI a combination of the differential measurement strategy with calibration and 
de-striping further reduce the effect. 



Table 2. Effect on Planck LFI maps of the main systematic effects 

30 GHz 44 GHz 70 GHz 

[MK] [mK] \p.K] 



Channel 


p-p 


rms 


P-P 


rms 


P-P 


rms 


1-Hz spikes . . 


4.00 


0.45 


1.51 


0.15 


2.56 


0.30 


Thermal fluct 














Back-end 


1.27 


0.11 


0.63 


0.05 


2.70 


0.24 


Front-end 


1.05 


0.23 


1.15 


0.22 


1.12 


0.21 


4K loads 


9.76 


0.98 


9.73 


0.98 


1.30 


0.16 


Total 


10.92 


1.10 


9.73 


0.98 


4.28 


0.45 



3. Component separation 

As widely discussed in the next section, the sky emission, at a given frequency, is 
a superposition of emission from various sources. The plausible contamination of 
the observable primary CMB by foreground emission has always been a source of 
concern for CMB observations. However, the level of foreground contamination, at 
high Galactic latitude and at frequencies between 50 and 200 GHz, is low enough 
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Fig. 1. Impact of main systematic effects in Planck LFI on temperature angular power spectrum 
(APS) as function of the multipole I. 

that for I less than about 2500 the temperature (or total intensity) APS, Cj T , 
can be accurately measured with only minor masking of the regions most contami- 
nated by foregrounds (Galactic interstellar medium (ISM) and bright extragalactic 
compact sources) . On smaller scales, emission from a background of blended faint 
extragalactic sources contributes a significant fraction of the observed powerS^HUl. 

For a sensitive mission such as Planck, foreground emission, rather than in- 
strumental noise, sets the limit of the accuracy of the measurement of the CMB 
APS. This limit depends on the effectiveness of any foreground-cleaning technique 
used to separate the primary CMB emission from foregrounds. Hence, the devel- 
opment, comparison, and optimization of component separation methods has been 

an important activity in the Planck Collaboration during the preparation of the 
missiojpMED. 

3.1. Modeling sky emission 

A multicomponent model of sky emission serves both as a framework for analyzing 
and interpreting the observations, and as a summary of our knowledge about the as- 
trophysical emitters. Recently, such a model of sky emission, the Planck Sky Model 
(PSM), has been put together for this purpose^. It is based on an underlying 
ACDM cosmological model with associated standard parameters^] In this frame- 
work, the emission of the CMB, of galaxy clusters, and of high-redshift galaxies in 
haloes of large-scale density contrast, is described on the basis of their statistical 
properties (angular power spectra of the CMB, cluster number counts, source num- 
ber counts and halo occupation distribution as a function of luminosity, of redshift, 
and of spectral energy distribution). Known bright radio and infrared sources are 
modeled using extrapolations of their measured fluxes at various frequencies. Galac- 
tic diffuse emission is modeled as a mixture of several components: synchrotron radi- 



d The spectral index and amplitude of scalar perturbation, n s and A s , the ratio between the 
amplitude of tensor and scalar perturbations, r = At / A s , the density parameters of matter, baryons, 
dark energy, f! m , Q^, Qa, the Hubble constant, Ho, etc., the default values of which are set to the 
current best fit. 
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ation from energetic electrons spiraling in the Galactic magnetic field, thermal dust 
emission, free-free (Bremsstrahlung) emission from the warm interstellar medium, 
the so-called 'anomalous' dust emission from small spinning dust grains. Molecular 
line emission could also contribute to the signal in some frequency bands, that are 
typically relatively wide in order to increase the measurement sensitivities. In par- 
ticular, the emission from the CO line is clearly seen in Planck data sets that allow 
to produce all-sky accurate maps of this signaP. 

3.2. Basics of component separation 

Consider a single pixel p in a set of maps observed at various frequencies, indexed 
by v. The signal observed at frequency v 1 in pixel p, can be written as 

x{v,P) = ^T j a%{v,p) 8i(p) + n(v,p), (1) 

i 

or, in vector-matrix format 

x(p)=A(p)s(p)+n(p). (2) 

If we know the frequency dependence A(p) of each foreground in pixel p, our problem 
is just to invert a (set of) linear system(s) to find the reference component templates 
s(p). In the limit where the instrumental noise n(p) is small, the inversion is imple- 
mented using the inverse of A (or, for non-square systems, the left pseudo-inverse 
Aj~f t = [A t A] _1 A t ). Otherwise, classical linear solutions such as least square (LS) or 
Wiener inversion can be used. They require, however, the prior knowledge of the 
covariance matrix R n (p) of the instrumental noise in each pixel, and possibly also 
of that of the signal, R s (p). The LS solution is 

? LS (p)=[A t R^ 1 A]- 1 A 1 R^ 1 x(p) ) (3) 

and the Wiener one is 

swiencr(p) = [A*R^A + R^" 1 ] -1 A*R^x(p), (4) 

where A* denotes the transpose of A. Most of the time however, neither the 'mixing 
matrix' A, nor the statistical properties of the signal, are known. Sometimes even 
the noise covariance is not well known. One must then find a way to estimate them 
from the data themselves prior to inverting the linear mixture and recovering an 
estimate of each of the components. 

3.3. Blind component separation 

Blind separation of linear mixtures is a classical field of research in signal and 
image processing. Typically, the data model is of the form of Eq. ^ , except that 
A does not depend on p, and noise is often not an issue. The main problem is 
then to decide how many components exist in the data, and determine A (or find a 
matrix W that inverts the system, without explicitly estimating A). The main idea 
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is to use the (assumed) statistical independence of the various components s. Blind 
component separation applied to CMB observations must address specific issues: 
error estimation, ill-conditioned covariance matrices, correlations between some of 
the components, spherical data sets. This has led the CMB community to adapt 
classical independent component analysis (ICA) methods for their analyses. 

The FastICA methocEOH aims at inverting the linear system using the matrix 
W that maximizes a measure of the non-Gaussianity of the component maps. The 
method, however, is not very effective at distinguishing Gaussian CMB from Gaus- 
sian noise. It also fails to exploit the strong spatial correlation of the CMB and of 
most of the diffuse foreground emission. 

Spectral matching ICA (SMICAJpMU, 

is a flexible method that maximizes the 
likelihood of a parametric model of A, R s and R s by minimizing the spectral mis- 
match between empirical and modeled second-order statistics of the observed maps. 
It is particularly useful for measuring a CMB APS, or parameters that model it, 
directly from multifrequency data. It has been used for predicting the errors on the 
tensor to scalar ratio r that can be reached by future CMB B-modes experimental. 
Very similar in spirit, although many implementation details vary, the correlated 
component analysis (CCA) method has been developed to deal specifically with 
correlated components^. Once second-order statistics of components and noise are 
obtained, they are used to invert the linear system. 

Neural networks provide another attractive solution for finding either coefficients 
that invert the linear mixture, or only those coefficients that recover the CMB 
specifically^. The method seems to perform well both on simulations^ and on 
real data^l. The impact of the training of the neural network, however, is hard to 
evaluate, and the propagation of errors not straightforward. These two limitations 
are serious for CMB data analysis, and may explain why the method has not received 
more attention so far. 

A completely different point of view is taken in |31| . Instead of learning (explic- 
itly or implicitly) the model (or part of it) in the data themselves, a parametric 
model of all relevant foregrounds is assumed a priori. Specifically, an amplitude 
and parameters defining an appropriate emission law are assigned to each emis- 
sion process in each sky pixel. The value of all parameters are then found using 
a Monte-Carlo Markov chain (MCMC) algorithm. This is clearly the appropriate 
approach for measuring efficiently parameters of a known model, and is hence of 
much interest for component separation. It is also very flexible, as the foreground 
model can be chosen freely. The main caveat is that while in theory this method 
provides a complete likelihood (and hence errors) for all parameters, in practice the 
main uncertainty is whether the assumed parametric model is correct. Goodness of 
fit is not a fully satisfactory criterion: as well known, it is always possible to fit a 
limited data set with a wrong model, provided the number of parameters used in 
the fit is large enough. The method is hence of interest only when many different 
channels of observation are available. An alternative implementation of this idea 
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has been proposed by [52] . 

Finally, yet another approach, based on a different optimization criterion to 
recover the mixing matrix A, has been proposed by [33] . It uses a likelihood penal- 
ization imposing a sparse representation of sky components in some over-complete 
dictionary of functions that serve as a redundant basis. The first implementation of 
the method, although promising and conceptually interesting, was not optimized for 
real-life data processing, due to lack of proper handling of different map resolutions 
and of the non-stationarity of the component emissions. This has been improved 
with a recent version that uses wavelet decompositions, allowing variations in both 
pixel and harmonic space of the linear combinations of the maps used for CMB 
recovery^. Such localization makes it possible to relax the restrictive condition 
that the emission of each component should be decomposable into the product of a 
pixel-independent emission law A{v) and a spatial template s{p). 



3.4. The internal linear combination and variants 

Bypassing in some way the need to rely on a specific model of all foreground emission 
(especially rigid linear mixtures in which A does not depend on p) is an appealing 
option. The frequency scaling of CMB anisotropics themselves being known to be 
the derivative with respect to temperature of a 2.725 K blackbody, independently 
of p, one may write a simplified model of the observed maps as 

x(p) = as(p) + n(p), (5) 

where s(p) is the CMB map, a the CMB frequency scaling, and all the unknown 
(or poorly known) noise and foreground contamination is dumped together into a 
single noise term n(p). The LS reconstruction of the CMB map is then 

This may seem impossible to implement without knowing R„ (which now includes 
unknown foregrounds correlated between channels). However, under the hypothesis 
that s(p) is not correlated with n(p), the covariance R x of the observed maps is R x = 
erjraa' + R„, where a\ is the variance of the CMB map. It is then straightforward 
to show, using the Woodburry inversion formula, that a'R^T 1 cx a*R^' 1 , and hence 
that the LS solution can be rewritten as 

a'R^xfo) 

SM = ^R?a- (?) 

This form is easily implemented using empirical estimates R^ of R^ obtained on 
the data themselves. This solution is also obtained as the constrained minimization 
problem of finding the linear combination w'x of the inputs that has minimum 
variance under the 'CMB preserving' condition w*a = 1. It is classically called the 
internal linear combination (ILC) method. 
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The ILC method, for its simplicity and robustness, has been used for the analysis 
of COBE-DMR, WMAP, and Planck datcHHEO. Variants that compute weights 
in different regions of pixel, harmonic or wavelet space for CMB temperature or 
polarization have been derivecP^S^. Extensions for recovering more components 
than just the CMB are discussed by [4T1 l42j . 

It is important to note that the ILC is prone to subtle biases, which must be 
understood and controlled for scientific analyses based on ILC maps. The first bias, 
a loss of some modes of the original CMB and hence of CMB power, is due to 
empirical correlations between the CMB and the contaminants, and is discussed at 
length in the appendix of [38] . The second is an amplification of calibration errors in 
the observed channels (or errors in the assumed frequency scaling of the component 
of interest), and is discussed in detail by [43] . 

3.5. Error assessment and masking 

One of the most crucial questions, once the component separation is performed, 
is the assessment of errors. How well does a method perform? While it is easy to 
propagate errors in a fit, the problem in component separation is that modeling 
errors dominate the uncertainties. Nonetheless, three approaches can give an idea 
of component separation performance. 

First, one can compare the results obtained with methods that are conceptually 
very different. If, however, results are very method-dependent, as usual, then one 
must either explain the differences and discard the method(s) thought to be less 
effective, or, for the post-analysis of the output map(s), mask (or flag as plausibly 
contaminated) any sky region in which agreement cannot be achieved. 

Second, one may test methods on simulations that are as realistic as possible. 
This has been one of the original motivations for the development of the PSM. 
However, the performance of some component separation methods is very sensitive 
on very subtle details about the sky emission. The refinement of the model and the 
separation of components are thus two complementary parts of a global, iterative, 
data analysis chain. 

Finally there is a third method, that permits to identify regions of the sky where 
the number of channels available is not sufficient to separate all emissions. Consider 
noisy observations of unspecified sky signals, x(p) = s(p) + n(p). Usually, the sky 
and noise components are pairwise de-correlated, and thus = R s + R n . We now 
suppose that R„ (instrumental noise only) is reasonably well known. Then we can 
whiten the observations (by multiplication by the square root of R n ). For the new 
data set, we have R^ = R s + Id. In the basis of diagonalisation of R x , the covariancc 
becomes R^ = A + Id. The number of (local) eigenvalues of R x significantly larger 
than unity is the dimension of the space spanned by measurable signal components. 
If all eigenvalues are larger than unity, then there are locally more independent 
astrophysical emissions than can be separated without external information. If, 
however, only some of the eigenvalues are significantly larger than unity, then in 
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Freq [GHz] 

Fig. 2. The expected spectral energy distribution (SED) for the sky averaged ZLE (blue line) 
and an hypothetical cold dust component with Td ust = 30 K, r = 1CT 7 ( green line) and their sum 
(red line) compared with the COBE/FIRAS spectrum (gray patch). 

principle the data set is redundant enough for blind component separation, for 
instance with methods such as SMICA and the ILC which, implemented in wavelet 
space, are expected to perform very satisfactorily. 

4. Diffuse foregrounds 

Except for the averaged (monopole) signal, microwave and sub-mm surveys are 
dominated at large angular scales by the diffuse signals from the Solar System and 
the Milky Way, emerging as foreground sources with patterns particularly prominent 
close to the ecliptic and Galactic plane, respectively, and, typically, remarkable up 
to few tens of degrees from them. The observer position plays a significant role in 
the study of the former (while it can be fully neglected for the latter, at least at the 
angular resolutions relevant here). For this reason, a special care in the application 
of the component separation methods described in previous section, if not a fully 
different approach, is required in this case, since, in general, the time dependence 
of the signal should be accounted for. 

4.1. Solar System diffuse emissions 

Solar System provides the most fore of all the foregrounds. In particular Zodiacal 
Light Emission (ZLE), i.e. the emission from Interplanetary Dust Particles (IDPs), 
dominated the sky signal at short wavelengths. At wavelengths shorter than 12 /im 
ZLE is mainly due to scattering of solar radiation, while at longer wavelengths 
thermal emission is the most important generation mechanisnEHIMl. ZLE is usually 
not accounted in CMB studies. In fact, since the ZLE flux below 1 THz decreases 
with i/ 4 (see Fig. [2| where the COBE/FIRAS spectrum is provided by 05]), i.e. as 
a modified blackbody with Tdust — 240 K, its contribution would be significantly 
smaller than the others foreground signals at CMB related frequencies. 

However Fig. [2] is based just on the observed emission of the IDPs at frequencies 
above 1 THz, which is dominated by the population of grains between Earth and 
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Jupiter orbits. Given dust orbiting the Sun is removed in short times by Pointing- 
Robertson decay, radiation pressure, grains-grains collisions and planetary encoun- 
ters a mechanism to continuously refurbish the population of IDPs is needed. In the 
inner Solar System the main contributors are comets and asteroids, even if some 
dust is supposed to come from interstellar space. Other sources of dust are know to 
be effective far from the Jupiter orbit, in particular the erosion and mutual collisions 
of trans-Neptunian objects (TNOs) and Kuiper Belt objects (KBOs) produced a 
second band of dust in the outer Solar System, detected by deep space probes. We 
will denote with KBOE the diffuse emission from this class of IDPs. 

Being at heliocentric distances much larger than IDPs responsible for standard 
ZLE, particles responsible for KBOE would be quite cold having Td us t ~ 30 — 60 K 
or lower. The optical depth below 1 THz is unknown. It depends on the balance 
between the mechanisms of dust production either mutual collisions between larger 
bodies or erosion by interstellar dust, the former tends to produce larger grains then 
the latter, the dust collision rate, the composition of dust, the geometrical distribu- 
tion of dust and the relative efficiency of production and destruction mechanisms. 
However it hardly will exceeds the optical depth for standard ZLE which is about 
10~ 7 . Two reasons make hard the detection of KBOE: (i) at high frequencies the 
ZLE emission would overwhelm the KBOE, as evident from the example of SED in 
Fig. [2j (n) the main method to separate the ZLE from background emission, basi- 
cally from Galactic dust, is to measure the seasonal dependence of ZLE signal for a 
fixed line-of-view, as the observer moves within the Solar System while surveying 
the sky. Such effect amounts to at most 7% of the emission, more than 90% of 
the modulation comes from dust within 3 AU from the Earth, so that the seasonal 
dependence from dust at heliocentric radii of 50 AU would be largely negligible. 

It is evident how Solar System diffuse emission may act as a foreground for CMB 
in two ways. The simplest case is when templates for Galactic and extragalactic 
emission at frequencies relevant for CMB observations are extrapolated from maps 
obtained at wavelengths where the ZLE is an important source of systematics. 
This is typical of CMB observations at frequencies higher than 100 GHz, where a 
template dust map is needed. Such maps are produced by extrapolating at CMB 
frequencies maps from observations made at several tens of //m where the ZLE is 
very strong and must be removed. In this case the ZLE will represent a sort of 
indirect foreground whose exact impact on CMB will depend. The other possibility 
is that of a direct impact on CMB missions of an un-removed component which is 
relevant at CMB frequencies, such as the KBOE. In this case the exact effect will 
depend on the level of contamination compared to the level of the CMB fluctuations, 
on the spatial pattern, and on its combination with the CMB. 

A simple argument allows us to estimate the effect of un-removed ZLE or KBOE 
on the CMB APS at different multipoles. Let us to consider a sky with just ZLE 
or KBOE and an observatory scanning the sky to form maps. Due to unavoidable 
geometrical constrains, the observed sky regions are more or less tightly correlated 
to specific positions in the Solar System, but it is not possible to grant a one- 
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to-one correlation, so that when timelines are co-added to form maps the result 
will be a distribution with small discontinuities, and in general different scanning 
strategies will create slightly different map;P^. However, if observations from many 
scans, spanning several years, are co-added, then a smooth map with a very strong 
planar symmetry about the ecliptic will be obtained. Denoting with Of jTnj ss the 
coefficients of the spherical harmonics (SH) expansion for such map, the planar 
symmetry assures that in a reference frame defined by this plane af i?Tij ss = for 
any odd I or for any m ^ 0, so that C| s = af q 2 / (2^ + 1) for even t and zeros for odd 
t. The combined map with CMB will be unchanged for odd I while will have in the 
same reference frame Ci = Cf MB + Cf s + a^ s Q af^ B / (21 + 1). So, depending on the 
signs of the ZLE or KBOE and CMB components, Ce can be smaller or larger than 
C™ B . Therefore, it can not be in principle excluded that part of the anomalies seen 
at low multipoles can be ascribed to some unknown and un-removed component of 
the ZLE or KBOB^ and/or to interplay between this foreground and un-removed 
dipole-like systematic effects^, especially at low I where a Solar System large scale 
diffuse emission should have the maximum power. 

4.2. Galactic emissions 

The wide frequency coverage of Planck when taken with relevant ancillary data 
provides a unique opportunity to characterize all the relevant Galactic foreground 
components. Of particular interest is the recently identified anomalous microwave 
emission (AME) due to spinning dust which has an important contribution at the 
lower Planck frequencies. Inclusion of this component has a domino effect on the 
spectrum of the other components, particularly at frequency v < 100GHz, where 
synchrotron and free-free emissions are particularly important. The emission from 
thermal (vibrational) dust dominates at v ^ 70GHz, but, although weak, it must 
be considered also at lower frequencies. It should be remembered that the minimum 
of the Galactic foreground to the CMB is in the range 60-100 GHz where each of 
these four components can have a small but significant contribution. 

4.2.1. Synchrotron emission 

Synchrotron emission originates in relativistic cosmic ray electrons spiraling in the 
Galactic magnetic field. The relativistic electrons are produced in the shocks asso- 
ciated with supernova explosions. The spectrum of the synchrotron radio emission 
is related to the energy spectrum of the relativistic electrons. Up to several GHz the 
brightness temperature spectral index is ~ —2.7^; above this frequency it steep- 
ens to —3.0 or more at the lower Planck frequencies^. Another characteristic of 
synchrotron emission is its linear polarization which is orthogonal to the magnetic 
field direction; this may be as high as 70 % for an aligned field with a brightness 
spectrum of —3.0. In the more tangled field environment on the Galactic plane the 
observed values are in the range 10-50%. 
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The current study of the emission of the plane in the inner Galaxy identified 
the synchrotron component by using component separation techniques. The low 
frequency data (0.4 to 2.3 GHz) revealed a narrow component in Galactic latitude 
with a FWHM of 1.6°. This newly identified component is also clearly identified in 
K and Ka band polarization data from WMAP; it has a brightness spectral index 
of —3.2^. This narrow distribution is the sum of the supernova remnants (SNRs) 
over the last 10 5 -10 6 years (the timescale of the SNRs before they expand into the 
broader latitude distribution. A similar latitude width is found for the normal (~1 
second period) pulsars; their ages are also ~ 10 5 -10 6 years. Over this timescale both 
the SNR shells and the pulsar proper motions will have taken them to a FWHP of 
1.6°, double the width of the nascent OB star distribution (0.9°). 

4.2.2. Free-free emission 

The free-free emission in the inner Galaxy arises from the ionized (electron) gas 
component produced principally by the hot O and B stars which are confined to a 
narrow latitude width of 0.9° FWHM. Cooler stars also contribute to the interstellar 
radiation field (ISRF) which is more diffuse. At intermediate and high latitudes the 
free-free emission is measured by the Ha spectral line. Even here a correction is 
needed to account for the absorption of Ha by dust. On the Galactic plane the 
dust obscuration is so great that the Ha emission line cannot be used. Here the 
radio recombination lines (RRLs) save the day. No dust absorption correction is 
required. However an electron temperature is needed to determine the emission 
measure (EM = n\V) in order to derive the corresponding continuum temperature 
at any frequency^. The brightness temperature spectral index is well determined at 
Planck frequencies; it is ~ —2.13 at 30 GHz. The electron temperature of the diffuse 
ionized gas appears to be similar to the average for the compact HII regiondHEH. 

The FWHM of the free- free (1.1°) is intermediate between that of OB stars 
(0.9°) and the neutral hydrogen (1.8°). This is not unexpected since the gas (HI, 
H2 and dust) density is greatest on the plane and also because the ionized emission 
is proportional to n^L. The free-free, along with the AME, dominates the emission 
on the plane in the inner Galaxy. 

4.2.3. Anomalous Microwave Emission (AME) 

AME is the recently identified emission component which is well-correlated with far- 
infrared (FIR) dust emission. It is produced by rapidly spinning small dust grains 
having an electric dipole moment- 4 - . Typical masses are ^50 atoms which in a dust 
cloud produce a spectrum which peaks in the range 15-50 GHz depending on the 
environment and radiation field. Planck has for the first time been able to define 
the shape of the spectrum on the high frequency side of the emission peak in a 
number of dust/molecular/HII regions, as shown in |55j . This work has provided a 
rich source of data to explore the emission mechanism in detail. 
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On the Galactic plane the AME spectrum can be estimated by applying com- 
ponent separation techniques to the strong signals measured here. The AME is the 
residual emission after the free-free, the synchrotron and the thermal dust have been 
accounted for. In the frequency range 20-40 GHz AME is comparable in brightness 
to the free-free for the inner Galactic plane from I = 300°-0°-60°. The latitude 
width of the emission at these frequencies is similar to that of the thermal dust. 

4.2.4. Thermal dust emission 

The FIR dust spectrum is due to the vibrational emission from dust grains heated 
by the ISRF. The peak in the emission is at a wavelength of ~60-100 microns. 
Averaged over the intermediate latitude sky the dust temperature is ~18K with a 
grey body slope in brightness of +1.7^. On the Galactic plane the dust temperature 
is somewhat higher at 20-24 K. The latitude width of the dust emission at say 
lOOmicrons is 1.2°, similar to that of CO (representing H2). The narrower width 
compared with HI is probably due to the higher dust temperature on the plane 
produced by the O and B stars. 

4.2.5. Emissions close to the Galactic plane 

We find a narrow latitude distribution on the Galactic plane for each of the four 
emission components, synchrotron, free-free, AME and thermal dust. Recent star 
formation over the last 10 5 -10 6 years in the dense gas regions on the plane is most 
likely the cause. 

Using precise full-sky observations from Planck, and applying several methods 
of component separation, the emission from the Galactic "haze" at microwave wave- 
lengths has been identified and characterized^. The haze is a distinct component 
of diffuse Galactic emission, roughly centered on the Galactic centre, and extends 
to |6| ~ 35° in Galactic latitude and \l\ ~ 15° in longitude. By combining WMAP 
and Planck data, [57] were able to determine the spectrum of this emission to high 
accuracy, unhindered by the large systematic biases present in previous analyses. 
The derived spectrum is consistent with power-law emission with a spectral index 
of —2.55 ±0.05, thus excluding free-free emission as the source and instead favoring 
hard-spectrum synchrotron radiation from an electron population with a spectrum 
(number density per energy) dN/dE ~ E~ 2 ' 1 . At Galactic latitudes \b\ < 30°, the 
microwave haze morphology is consistent with that of the Fermi gamma-ray " haze" 
or "bubbles" (see also |58J, indicating that we have a multi- wavelength view of 
a distinct component of our Galaxy. Given both the very hard spectrum and the 
extended nature of the emission, it is highly unlikely that the haze electrons result 
from supernova shocks in the Galactic disk. Instead, a new mechanism for cosmic- 
ray acceleration in the centre of our Galaxy is implied. 

The wide frequency coverage of Planck, which includes polarization, allows the 
spectrum of each component to be determined unambiguously. Polarization data 
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from Planck are awaited with considerable interest. 

5. Extragalactic radio and far IR sources at mm/sub-mm 
wavelengths 

The Planck Early Release Compact Source Catalogue (ERCSC)P^ - the first com- 
plete full-sky catalogue of bright sub-millimeter extragalactic compact sources - pro- 
vides positions and flux densities of hundreds of "radio" sources (intermediate to 
high-redshift Active Galactic Nuclei (AGN)) and of thousands of "far-IR" sources 
(low-redshift dusty galaxies) detected in each of the nine Planck frequency maps 
during the first 1.6 Planck full-sky surveys. As shown in [55], their Table 1, the 
full-sky surveys of the Planck satellite are - and will be, for years to come - unique 
in the millimeter, at A < 3 mm, and sub-millimeter domains. Thanks to this huge 
amount of new data it is thus possible to investigate the SEDs of extragalactic 
point sources in a spectral domain very poorly explored before and, at the same 
time, their cosmological evolution, at least for some relevant source populations. 

5.1. Radio sources: "blazars" 

The most recent estimates on source number counts of extragalactic radio (syn- 
chrotron) sources up to ~ 50 — 70 GHz, and the optical identifications of the 
corresponding point sources (see e.g. [60]), show that these counts are dominated 
by radio sources whose average spectral index is "flat", i.e., a ~ 0.0 (with the usual 
convention S v oc v a ). This result confirms that the underlying source population 
is essentially made of Flat Spectrum Radio Quasars (FSRQ) and BL Lac objects, 
collectively called "blazars"^] with minor contributions coming from other source 
population jj 62 ll 63 l . At frequencies > 100 GHz, however, there is now new information 
for sources with flux densities below about 1 Jy, coming from the South Pole Tele- 
scope (SPT) collaboration's, with surveys over 87 deg 2 at 150 and 220 GHz, and 
from the Atacama Cosmology Telescope (ACT) survey over 455 deg 2 at 148 GHz^. 

To study the spectral properties of the extragalactic radio sources in the Planck 
ERCSd^Sl used a reference 30 GHz sample above an estimated completeness limit 
Sum — 1-0 Jy. In this sample, the 30-143 GHz median spectral index is in very 
good agreement with the one found by Marriage et alP^ for their bright (S^ > 50 
mJy) 148 GHz-selected sample with complete cross-identifications from the Aus- 
tralia Telescope 20 GHz survey, i.e a^f — —0.39 ± 0.04. In the whole, the results 
of show that in their sample selected at 30 GHz a moderate steepening of the 
spectral indices of the radio sources at high radio frequencies, i.e. > 70 — 100 GHz, 
is clearly apparent. It has also been shown by |66j that differential number counts at 
30, 44, and 70 GHz are in good agreement with those derived from WMAP dataP^ 



c Blazars are jet-dominated extragalactic objects characterized by a strongly variable and polarized 
emission of the non-thermal radiation, from low radio energies up to high energy gamma rays; see 
e.g. [61]. 
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at nearby frequencies. The model proposed by de Zotti et alP^in 2005 is consistent 
with the present counts at frequencies up to 70 GHz, but over-predicts the counts 
at higher frequencies by a factor of about 2.0 at 143 GHz and about 2.6 at 217 GHz. 
As reminded before, the analysis of the spectral index distribution over different 
frequency intervals, within the uniquely broad range covered by Planck in the mm 
and sub-mm domain, has highlighted an average steepening of source spectra above 
about 70 GHz. This steepening accounts for the discrepancy between the model 
predictions of de Zotti et alP^ and the observed differential number counts at HFI 
frequencies. 

Recently, a successful explanation of the change detected in the spectral behav- 
ior of extragalactic radio sources (ERS) at frequencies above 70-80 GHz has been 
proposed by Tucci et alP^. By applying the KonigP^ model for the emission in the 
inner jets of blazars, |68) makes a first attempt at constraining the most relevant 
physical parameters that characterize the emission of blazar sources by using the 
number counts and the spectral properties of extragalactic radio sources estimated 
from high-frequency radio survey^] As noted before, a relevant steepening in blazar 
spectra with emerging spectral indices in the interval between —0.5 and —1.2, is com- 
monly observed at mm/sub-mm wavelengths. Tucci et al.^ interpreted this spectral 
behavior as caused, at least partially, by the transition from the optically-thick to 
the optically-thin regime in the observed synchrotron emission of AGN jetsP, giv- 
ing rise to a "break" frequency, vm, typically in the range between 50-2000 GHz, at 
which the synchrotron spectrum of jets bends dowrjs] On the whole, the results of 
[B"Bj imply that the parameter tm should be of parsec-scales, at least for FSRQs, in 
agreement with theoretical predictions^, whereas values of r^/ <C lpc should be 
only typical of BLLac objects or of rare, and compact, quasar sources. 

5.2. Far IR sources: local dusty galaxies 

The analysis done by |72) presented the first results on the properties of nearby 
galaxies using ERCSC data. From reliable associations between Planck and IRAS, 
they selected a subset of 468 for SED studies, namely those with strong detections in 
the three highest frequency Planck bands and no evidence of cirrus contamination. 
This selection has thus provided a first Planck sample of local, i.e. at redshift < 0.1, 
dusty galaxies R The analysis of SEDs of these local galaxies^! has confirmed the 



The main goal of 68] was to present physically grounded models to extrapolate the number counts 
of ERS, observationally determined over very large flux density intervals at cm wavelengths down 
to mm wavelengths, where experiments aimed at accurately measuring CMB anisotropics are 
carried out. 

s Bascd on published models, Tucci et al.EH estimated the value of the frequency um (and of 
the corresponding radius r^) at which the break occurs on the basis of the ERS flux densities 
measured at 5 GHz and of the most typical values for the relevant physical parameters of AGN 
jets. 

h This sample is very important for determining their emission properties and, in particular, the 
presence of different dust components contributing to their sub-mm SEDs. 
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presence of cold dust in local giant and, largely, in dwarf galaxie^] In [72] it is also 
found that some local galaxies are both luminous and cool, with properties similar to 
those of the distant SMGs uncovered in deep sub-mm surveys. The main conclusion 
of [72] is that cold (T < 20 K) dust is thus a significant and largely unexplored 
component of many nearby galaxies and that there is a new population of cool 
sub-mm galaxies, showing the presence of even cooler dust grains, with estimated 
temperatures of ~10-13 K. 

Very recently, using selected samples from the first Planck 1.6 full-sky surveys, 
i.e. the Planck ERCSC, [73] derived number counts of extragalactic point sources 
from 100 to 857 GHz (3 mm to 350 fim). More specifically, for the first time, num- 
ber counts have been provided of synchrotron dominated sources (blazars) at high 
Planck frequencies (353 to 857 GHz) and of dusty galaxies at lower frequencies (217 
and 353 GHz). Planck number counts are found to be in the Euclidean regime in 
this frequency range, since the ERCSC comprises only bright sources (S > 0.3 Jy). 
The estimated number counts appear generally in agreement with other data sets, 
when available (see [73] for more details). 

These new estimates of number counts of synchrotron and of dust-dominated 
extragalactic sources allowed new constraints to be placed on cosmological evolution 
models which extend their predictions to bright flux densities, i.e. S > 1 Jy. A very 
relevant result is that the most successful model of Tucci et alP^ is performing par- 
ticularly well at reproducing the number counts of synchrotron-dominated sources 
up to 545 GHz. On the contrary, |73j highlights the failure of many models for num- 
ber count predictions of dusty galaxies to reproduce all the high-frequency counts. 
The likely origin of these discrepancies is an inaccurate description of the galaxy 
SEDs used at low redshift in these models. Indeed a cold dust component, detected 
by [72] , is rarely included in the models of galaxy SEDs at low redshift. On the 
whole, these results already obtained by the exploitation of the Planck ERCSC data 
are providing valuable information about the ubiquity of cold dust in the local Uni- 
verse, at least in statistical terms, and are guiding to a better understanding of the 
cosmological evolution of extragalactic point sources at mm/sub-mm wavelengths. 

5.3. Nearby galaxies: the case of M31 

As discussed above, WMAP and Planck data can be used to get information about 
point like astrophysical sources (see e.g. |74j and references therein), including 
nearby galaxies. Recently, the 7-year WMAP data have been used to trace the 
disk and the halo of the M31 galaxy-^. Unexpectedly, an asymmetry in the mean 
microwave temperature towards both the M31 disk and halo, along the direction of 



'The SEDs are fitted using parametric dust models to determine the range of dust temperatures 
and emissivities. They found evidence for colder dust than has previously been found in external 
galaxies, with temperatures T < 20 K. Such cold temperatures are found by using both the 
standard single temperature dust model with variable emissivity /?, or a two dust temperature 
model with p fixed at 2. 
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the M31 rotation, has been found. The maximum temperature contrast (see Fig. [3] 
reprinted from [75]) is about 130 /iK/pixel (or about 200 /jK/ 'pixel if the M31 Bulge 
is excluded). This temperature asymmetry, similar in the three WMAP bands W, 
V and Q, is very likely induced by the Doppler shift effect due to the M31 disk 
rotation speed. A similar effect is clearly visible also towards the M31 halo up to 
about 120 kpc from the M31 center with a peak value of about 40 /iK/pixel. 

The robustness of this result has been tested, both for the M31 disk and halo, by 
considering 500 randomly distributed control fields in the three WMAP bands and 
also by simulating 500 sky maps (see [76] for more details) . CMB maps are simulated 
by assuming AT(n) = ATcMB(n) <8> B(n) + N(n), where ATcmb is a realization 
of the Gaussian CMB field, N(n) is the pixel noise and B(n) is the proper beam 
of the experiment. Using the synfast routine of HEALPix-^ with the best-fit power 
spectrum constrained with BAO and Hq, as given by the WMAP Collaboration, 500 
realizations of the CMB sky were made. The maps have been then convolved with 
the WMAP beams for W, V, and Q bands, respectively, taking into account the 
convolution with the beam function of the experiment and randomly extracting the 
noise value from a normal distribution with a = a^j \J N i, s . The statistical analysis 
shows that there is a probability below about 1% that the temperature asymmetries 
both in the M31 disk and halo are due to random fluctuations of the CMB signal. 

The degree to which galactic halos rotate with respect to the disks is a difficult 
task to be investigated. In this respect, the methodology of using CMB data to 
probe both the disk and the halo of M31, even if with the limitation of the presently 
available data, may suggest a novel way of approaching this problem especially in 
view of the high accuracy CMB measurements with the Planck satellite. 




Fig. 3. The detailed geometry (up to 8°) used in the analysis is shown. The different pixel colors 
indicate the difference of the CMB temperature with respect to the average temperature set to 
zero. Red means positive excess and goes up to a maximum of 150/xK while blue means lower 
temperature and goes up to — 125/iK. The left image shows the real WMAP W band map, while 
the right image shows the geometry used in the analysis and the average temperature in the two 
sides of the M31 disk and halo in false colors. It shows in a single glance that one side of both the 
M31 disk and halo is hotter with respect to the other side. 
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5.4. Cosmic Infrared Background anisotropics 

The Cosmic Infrared Background (CIB) is the relic emission, at wavelengths larger 
than a few microns, of the formation and evolution of the galaxies of all types, 
including Active Galactic Nuclei (AGN) and star-forming systemJ^^j The CIB 
accounts for roughly half of the total energy in the optical/infrared Extragalactic 
Background Light (EBL)P^, although with some uncertainty, and its SED peaks 
near 150 /im. Since local galaxies give rise to an integrated infrared output that 
amounts to only about a third of the optical one^l, there must have been a strong 
evolution of galaxy properties towards enhanced far-IR output in the past. There- 
fore, the CIB, made up by high density, faint and distant galaxie^Jis barely resolved 
into its constituents. Indeed, less than 10% of the CIB is resolved by Spitzer at 160 
am^, - 10% by Herschel at 350 /xm^ and - 16% by the SCUBA-2 Cosmology 
Legacy Survey (S2CLS) at 450 /im^l. With the advent of large area far-IR to mil- 
limeter surveys (Herschel, Planck, SPT, and ACT), CIB anisotropics thus constitute 
a new tool for structure formation and evolution studies. 

Because the clustering of dark matter is reasonably well understood, observations 
of anisotropies in the CIB constrain the relationship between dusty, star-forming 
galaxies at high redshift, i.e. z > 2, and the underlying dark matter distribution. 
The APS of CIB anisotropies has two contributions: a white-noise component caused 
by shot noise and an additional component caused by spatial correlations between 
the sources of the CIB. Correlated CIB anisotropies have already been measured by 
many space-borne as well as ground-based experiments (see [87] for more details). 
On small angular scales (£ > 2000), they measure the clustering within a single 
dark matter halo and, accordingly, the physics governing how dusty, star-forming 
galaxies form within a halo. On larger angular scales, i.e. 200 < £ < 2000, CIB 
anisotropies measure clustering between galaxies in different dark matter halos. 
These measurements primarily constrain the large-scale, linear bias, b, of dusty 
galaxies, which is usually assumed to be scale-independent over the relevant range. 

Thanks to the exceptional quality of the Planck data, |87j were able to measure 
the clustering of dusty, star-forming galaxies at 217, 353, 545, and 857 GHz with 
unprecedented precision. After careful cleaning, based on suitable templates and 
Planck maps, they obtained CIB anisotropy maps that reveal structures produced 
by the cumulative emission of high-redshift, dusty, star-forming galaxies. The power 
spectra of the latter maps were then computed with high signal-to-noise ratio over 
the range 200 < I < 2000 by [87]. These measurements compare very well with 



J An important goal of studies about galaxy formation has thus been the characterization of the 
statistical behavior of galaxies responsible for the CIB - such as the number counts, redshift 
distribution, mean SED, luminosity function, clustering - and their physical properties, such as 
the roles of star-forming vs. accreting systems, the density of star formation, and the number 
density of very hot stars. 

k The CIB records much of the radiant energy released by processes of structure formation occurred 
since the decoupling of matter and radiation, four hundred thousand years after the Big Bang, 
when the CMB was produced. 
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previous measurements at higher ^] Moreover, from Planck data alone [57J could 
exclude a model where galaxies trace the (linear theory) matter power spectrum 
with a scale-independent bias: that model requires an unrealistic high level of shot 
noise to match the small-scale power they observed. Consequently, an alternative 
model that couples the dusty galaxy, parametric evolution model of [55J with a halo 
model approach has been developed (see, again, [57] for more details). Characterized 
by only two parameters, this model provides an excellent fit to the measured CIB 
anisotropy APS for each frequency treated independently. 

6. Clusters of galaxies and their cosmological implications 

The observation of clusters of galaxies through the Sunyaev-Zel'dovich (SZ) effect, 
the inverse Compton scattering of cosmic microwave photon by hot intra-cluster 
electrons®, have proven to be an efficient way to search for new clusters^H2ED. 

The Planck satellite has been observing clusters of galaxies via the measurement 
of the SZ effect over the whole sky since August 2009. Although, its spatial resolution 
is moderate with respect to ground based SZ surveys (see e.g. [55] [HS]), it possesses 
a unique nine-band coverage from 30 to 857 GHz and, most crucially, it covers an 
exceptionally large survey volume. Indeed Planck is the first all-sky survey capable 
of blind cluster detections since the ROSAT All-Sky Survey (RASS, in the X-ray 
domain). Early Planck results on galaxy clusters were recently published in [911 192[ 
|9"51 [Ml 1931 191)] . These results include the publication of the high signal-to-noise 
ratio (S/N > 6) Early SZ (ESZ) cluster sampleP. 

6.1. Planck SZ clusters 

Using this specific SZ signature, Planck was designed to be able to detect numerous 
clusters®. Unfortunately, not all are showing up as Abell 2319. The signal is indeed 
quite weak and is contaminated by foregrounds (our Galaxy, and nearby radio/IR 
galaxies) and backgrounds (CMB and CIB). As described later, the published Planck 
clusters have a signal-to-noise ratio (S/N) greater than 6. This means that the S/N 
per frequency is of the order of 1. This has lead us to develop a specific approach 
for detecting, validating and confirming clusters. 

We use a multi-matched filter (MMF) method® to detect the clusters. It is tak- 
ing advantage of the spectral signature (SZ signature without relativistic effects) 
and the spatial signature (universal spherical profile from X-ray REXCESS obser- 
vations®) of the clusters detected by Planck. As optimal as the method can be, 
a process of validation is still necessary to remove false detections. This is done 
in two steps. First a cross-check with internal Planck catalogues (cold cores, solar 
system objects, bad pixels) is performed, then cross-checks with existing external 



'The SED of CIB anisotropics is not different from the CIB mean SED, even at 217 GHz. This is 
expected from the model of | 88 | and reflects the fact that the CIB intensity and anisotropies are 
produced by the same population of sources. 
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Fig. 4. Distribution on the sky of the Planck ESZ clusters (the signal has been amplified to be 
seen). 

catalogues and data (SDSS, RASS) are performed to classify the known clusters and 
the new candidate clusters. Finally, follow-up observations has been done in optical, 
SZ and mainly in X-ray with XMM-Newton, to confirm our candidate clusters. 

6.1.1. Planck Early SZ cluster sample 

These detection, validation, and confirmation steps have lead to the production of 
the Planck Early SZ Cluster sample (ESZ). It contains 199 clusters, 10 of which, 
confirmed by XMM-Newton validation progra m 91 !! 100 ! 101 !, have a S/N <6. The 189 
clusters with S/N greater than 6 are divided in 169 known clusters (in X-ray, optical 
or SZ) and 20 new Planck clusters. At the time of the release only 11 were confirmed 
by XMM-Newton. Since then, 6 more have been confirmed by SPT and AMlIMIHl. 
The sample is availably™] as part of the Planck Early Release Compact Source 
Catalogue (ERCSC)P^. The distribution on the sky of these clusters is shown in 
Fig. [4] (reprinted from |104j ). 

The ESZ clusters have relatively low redshift; 86% of them have z < 0.3. Their 
masses span more than a decade up to 1.5 • 10 M so i, and a large fraction of new 
Planck detected clusters are massive (> 9 • 10 14 M soi ). Planck has thus a unique 
capability to detect the rarest and most massive clusters over the full sky. 

6.1.2. SZ clusters properties 

Observing galaxy clusters in SZ opens a new observational window to understand 
the clusters themselves and the evolution of our Universe. Planck has detected new 
clusters, sometimes massive. Why have they not been detected already in X-ray? 
Is this a new population of clusters, or the gas (responsible for both X-ray and SZ 
emissions) properties differ from what we think? As massive objects, clusters are 



rssd.esa.int /Planck 
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Fig. 5. Left panel: Planck SZ and XMM X-ray images of PLCKG214.6+37.0. Right panel: 
electronic density profiles of Planck and REXCESS clusters. 

sensitive to cosmological initial conditions and cosmic evolution. To use clusters for 
cosmological studies we need to relate their mass to our observation (SZ effect or 
Y-parameter). But is SZ effect a good proxy for the mass? How does the SZ signal 
relates to the X-ray luminosity, to the richness of clusters? The Planck ESZ clusters 
and Planck data are and will help in answering these questions. 

6.1.3. New Planck clusters 

The new Planck confirmed clusters have been compared with REXCESS X-ray de- 
tected clusters. Planck clusters show a more complex morphology, being sometimes 
really diffuse, extended, disturbed, and also double or triple. For the same given 
mass, they are also sub-luminous in X-ray compared to the REXCESS ones. Their 
electronic density profiles is on average lower in the center than the REXCESS ones 
(see Fig.[5j reprinted from [HI]) - Multi-wavelenght studies will help understand these 
properties. For example, |105] have observed one the XMM confirmed Planck new 
clusters and found radio arcs. Such findings, revealing shocks and/or merger, would 
imply higher temperature areas, that could enhance the SZ signal and explain why 
these clusters are seen in SZ and not in X-ray. More dedicated multi-wavclcnght 
studies are thus needed to better understand these clusters. 

6.2. Baryons in clusters of galaxies as seen in the Planck survey 

The total SZ signal is closely related to the cluster mass (see e.g. |106j ). and its 
surface brightness insensitive to distance. Therefore, SZ surveys can potentially 
be used to built unbiased close to mass selected cluster samples up to high red- 
shift. These scaling relations also bear the imprint of all gravitational and non- 
gravitational physical processes at play in the process of structure formation and 
evolution. Therefore such SZ samples of galaxy clusters will be of tremendous help 
for structure formation studies and to provide CMB independent cosmological con- 
straints (see e.g. [1071 11081 1109] . However, this requires a precise understanding of 
the statistical properties of the cluster population and furthermore a precise cali- 
bration of scaling relations between clusters physical properties and their mass. In 
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the following, we focus on the current results on SZ scaling relations with respect to 
PlancKs results. From three different approaches, we have brought tight constraints 
on the scaling relations between the SZ signal and clusters physical quantities. 

The statistical combination of ~ 1600 MCXC clusters at 0.01 < z < lEfl] with 
the all-sky Planck data led to a precise measurement of the correlation between the 
SZ signal and the X-ray luminosity. Averaging SZ fluxes in bins of X-ray luminosity, 
Lx, we detected the SZ signal at very high significance. This Planck observed signal 
is consistent with X-ray based predictions over two decades in X-ray luminosity, 
down to Lx — 10 43 erg/s <J L 500 E(z)~ 7 / 3 ^ 2 x 10 45 erg/s. We found no deficiency 
in SZ flux with respect to the X-rays within i?soo- This results underlines the 
robustness and consistency of our overall view of intra-cluster medium properties 
(left panel of Fig. |6j reprinted from Fig. 4 in :93J). This analysis fully agrees with 
the similar study carried on beforehand on the WMAP-5 data by [19] . 

Moreover, it is also consistent with the more in-depth investigation of the local 
scaling relations conducted over a sample of 62 massive known clusters detected by 
Planck at a high signal-to- noise ratio and with archival XMM-Newton dataP^. This 
analysis has allowed us to investigate the scaling relations between the SZ signal, 
D \ 5^00) an d the X-ray-derived properties (i.e., gas mass -Mg 5 oo, temperature Tx, 
luminosity i500jo.i-2.4i kcVi SZ signal analogue Yx.500 — A'ig.500 x Tx, and total mass 
-^500 )■ The derived results are in excellent agreement with both X-ray predictions 
and recently-published ground-based data derived from smaller samples (middle 
panel of Fig. |6j reprinted from Fig. 4, left panel, in [M]; see [11 H 1112] ). 

Finally, as stressed in the previous section, the new clusters detected by Planck 
follow the same scaling relations out to z = 1 without significant deviation from 
self-similar evolution, exhibiting an equivalent agreement between their SZ and X- 
ray properties as show on the right panel of Fig. [6] (reprinted from Fig. 7 in |101j : see 
[9"31 UOOj ). This behavior is seen down to an SZ signal of Y500 ~ 3 x 10~ 4 arcmin 2 . 
Below this threshold, we reach the current detection limit of Planck and Malmquist 
bias clearly appears (for details see |101j ). 

As pointed out recently by |113) . some observational and/or survey biases may 
arise from the combination of different effects and systematic biases. On the side 
of observable biased, one can mentioned the well known hydrostatic equilibrium 
hypothesis which biases the X-ray masses low with respect to the true mass by 
10-20%Elliii5l. The richness and weak lensing mass estimators from optical ob- 
servations require a better control of their individual and statistical measure- 
ments^illliiiltliHI . Finally for SZ measurements, cross-calibration between Planck, 
SPT and ACT measurement are certainly needed to further lower the photometric 
uncertainties and assess possible SZ flux measurement systematics. 

Some of the aforementioned biases have been investigated on an individual clus- 
ter basis in |120j . where the relation between the Planck SZ signal and the mass 
was studied using total masses derived from both weak lensing (WL) measure- 
ments r Hilli2 2 J and from Xray data assuming hydrostatic equilibrium (HE; XMM- 
Newton observations). While the A/wl — relation was consistent with previ- 
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Fig. 6. Left panel: Scaling relation between Planck SZ measurements and X-ray luminosity for ~ 
1600 MCXC clusters. Individual measurements are shown by the black dots and the corresponding 
bin averaged values by the red diamonds. Thick bars give the statistical errors, while the thin bars 
are bootstrap uncertainties. The X-ray based model is shown as a solid blue line, and the bin- 
averaged SZ cluster signal it predicts is shown by the blue stars. The red dot-dashed line shows the 
best fitting power-law to the data. Middle panel: SZ flux vs prediction from X-rays. Blue stars 
indicate cool core systems. The dashed line is the prediction from REXCESS X-ray observations. 
Right panel: Relation between apparent SZ signal (Y500) and the corresponding normalized Yx 
parameter. Black points show clusters in the Planck-ESZ sample with XMM-Newton archival data; 
green and red points represent Planck clusters confirmed with XMM-Newton validation program. 
The blue lines denote the Y500 scaling relations predicted from the REXCESS X-ray observations. 
The grey area corresponds to median Y500 values in Yx bins with ±1ct standard deviation. 



ous measurements using WL masses^, there was an offset in normalization with 
respect to the relation obtained using HE X-ray mass measurements. Since both 
the SZ measurements and the HE X-ray masses were consistent with our previ- 
ous work, we concluded that the normalization offset in the Mwl — is due 
to the X-ray masses being ~ 20 per cent higher than the WL masses. This is an 
unexpected result, given that simulations generally predict that HE X-ray masses 
should be smaller than WL masses owing to a the neglect of pressure support from 
bulk gas motions in the HE mass equation. Further investigation showed that the 
discrepancy is enhanced in dynamically disturbed systems and appears correlated 
with differences in mass concentration and the offset between the X-ray peak and 
the BCG position (the centers used for the X-ray and Wl mass determinations, 
respectively). More work is clearly needed, as discussed extensively in |120j . These 
remaining 10-20% inconstancies in scaling relations between SZ, X-ray and optical 
data are at hand. They need to be further investigated and quantified in order to 
reach an holistic view of the galaxy cluster properties. 

Together with the SPT and ACT telescopes, the Planck survey is shading new 
light on the population of galaxy clusters complementing our existing view of the 
ICM hot gas from the X-ray observations with high precision multi-frequency sub- 
millimeter to centimeter measurements. The scaling properties of the SZ signal 
together with other cluster observables have been investigated with various means 
and methods. Well constrained scaling relations between the SZ and X-ray mea- 
surements have been derived, with high precision calibration for the Y500 — ix,5oo> 
^500 — ^x,500 an d Y500 — M relations. These results emphasize the well consistent 



15, 2013 2:10 WSPC/INSTRUCTION FILE CBsessMG13review 



26 C. Burigana, R. D. Davies, P. de Bernardis, et al. 

picture we have of the ICM at least within i?5oo- Further Planck intermediate results 
are currently being published providing further insights on the clusters of galaxies. 

6.3. Neutrino mass from SZ surveys 

We explore here the possibility of setting useful constraints on the total neutrino 
mass from cluster number counts obtained by the ongoing Planck/SZ and future 
cosmic-variance-limited surveys. The precision with which this mass can be deter- 
mined from SZ number counts is limited mostly by uncertainties in the cluster 
mass function and intra-cluster gas evolution. We find that projected Planck/SZ 
cluster counts could yield the total neutrino mass with a (la) uncertainty of 0.06 
eV, assuming the mass is in the range 0.1 — 0.3 eV, and the survey detection limit 
is set at the 5a significance level. Based on expected results from future cosmic- 
variance-limited (CVL) SZ survey, we predict a la uncertainty of 0.04 eV, a level 
comparable to that expected when CMB lensing extraction is carried out with the 
same experiment. 

CMB measurements already placed meaningful upper limits on the total neu- 
trino mass from its impact on the early integrated Sachs- Wolfe effect. The energy 
scale of recombination, ~ 0.3 eV, sets this upper limit; if the total neutrino mass 
is larger than this value, then neutrinos are non-relativistic and do not contribute 
to the decay of gravitational potentials shortly after recombination. If, on the other 
hand, the total mass is lower they constitute a relativistic component that con- 
tributes to the decay of linear gravitational potentials, changing the temperature of 
the CMB towards these gravitational wells. 

Applying optimal estimators to CMB temperature and polarization maps one 
can recover the lensing potential to the precision that will allow constraining the 
total neutrino mass to the 0.04 eV level (see [123 ) with a cosmic-variance-limitcd 
(CVL) CMB experiment, assuming full-sky coverage, no foregrounds, and no source 
of non-Gaussianity other than the lensing of the CMB. In practice, it is unlikely that 
all these conditions will be fully satisfied and in that sense the frequently-quoted 
value 0.04 eV is likely to be unrealistic. 

Cluster number counts are yet another useful probe of neutrino masses. This is 
due to the fact that typical cluster scales are much smaller than the ~ 150 Mpc 
scale of linear dark matter halos that lens the CMB. In addition, cluster number 
counts are exponentially sensitive to a(M, z), the rms mass fluctuation on a cluster 
mass scale M at redshift z, and since a(M, z) itself is exponentially sensitive to 
neutrino mass (via the growth function), this implies that cluster number counts 
should be a rather sensitive probe of neutrino masses (see e.g. |1241 1125| L 

We further explored the ability to strengthen the constraints on the neutrino 
mass from cluster number counts, and extended our forecast to additional surveys. 
This was done by parameterizing uncertainties in the halo mass function, which 
is the dominant source of modeling uncertainties. The shape and normalization 
of the mass function reflect the details of the growth of density fluctuations, and 
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the nonlinear collapse and merger of sub-structures, whose hierarchical evolution 
can be best studied by state-of-the-art, large-volume hydrodynamical cosmological 
simulations. Currently available numerical codes predict a range of mass functions; 
this indeterminacy largely sets the precision limit of forecasting the total neutrino 
mass from cluster SZ number counts and power spectra. Additionally, we have 
accounted for cluster sample variance errors (in addition to Poissonian noise), a 
more realistic intra-cluster gas profile, as well as gas evolution with cluster mass 
and redshift. A full description of this work and results is given by J126j . 

Our analysis shows that from cluster number counts alone (and priors based on 
measurements of the primary CMB APS and the HST prior on Hq), the uncertainty 
in the determination of the total neutrino mass can be limited to the ~ 0.04 — 
0.06 eV range, depending on the details of the SZ cluster surveys and the fiducial 
neutrino mass. CMB anisotropy data combined with Planck cluster number counts 
are predicted to reach a level uncertainty at the higher end of this interval, whereas 
a CVL SZ survey is predicted to yield the somewhat higher precision corresponding 
to the lower end of this mass uncertainty interval. These results are based on the 
mass function by |127j , whose parameter values were assumed to have uncertainties 
that are higher by 10% than those specified there. 

7. Selected topics in CMB studies 

The release of first cosmological products and papers from the Planck mission, 
waited for early 2013, will have a strong impact for cosmology in the coming decades, 
following the very important results from WMAP and recent ground-based projects 
together covering a wide multipole range. In this section we discuss three very dif- 
ferent topics: the first is connected to fundamental physics results expected in next 
times from the Planck mission; the second, regarding the polarization imprints in- 
duced by galaxy clusters and filaments, is relevant in particular for future high 
resolution ground-based experiments; the last concerns the information on primor- 
dial power spectrum at extremely high wavenumbers that could be derived from 
next generations of CMB spectrum missions, thus linking absolute measures of the 
CMB monopole, i.e. the largest angular scale, to small scale phenomena. 

7.1. Fundamental physics from CMB Parity analyses 

The statistical properties of the CMB pattern may be used to constrain Parity 
(P) symmetry. Parity violations arise in several models: as modification of elec- 
tromagnetisnP^Hi22l or as modification of the standard picture of the Inflationary 
mechanism, where P is broken due to primordial (chiral) gravitational waves^UtlM'. 
Both of these scenarios predict non null cross-correlations between gradient and curl 
modes and scalar and curl modes in the CMB polarization pattern. However, chiral 
gravity induces such correlations at the CMB last scattering surface whereas cos- 
mological birefringence induces them by rotating the polarization plane during the 
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CMB photon journey from its last scattering to ui^il. We focus here mainly on cos- 
mic birefringence, reporting findings from Gruppuso et al. 2011'^'. In addition, we 
discuss the claimed P anomaly found at large angular scales in the anisotropy inten- 
sity spectrum of the WMAP data, first claimed by Kim and Naselsky in 201CP2SHH11 
The latter is dubbed a parity anomaly in view of an observed discrepancy (in power) 
among even and odd multipoles, which behave differently under P transformation. 
However, there is no sound theoretical framework that could explain such a mis- 
match. If the effect is indeed due to fundamental physics, its appearance at large 
angular scales naturally suggests the possibility that a P violating mechanism is 
involved during an early phase of the universe. Other explanations exist: for a more 
conservative approach see [137] where it is conjectured that we may live in a spe- 
cial location of the universe, such that translational invariance is violated at scales 
larger than ~ 4 Gpc. 



7.1.1. TT Parity anomaly 

All-sky temperature maps, T(n), are usually expanded in terms of spherical harmon- 
ics Yt m (n), with n being a unit vector or direction on the sky, completely specified 
by a couple of angles (9,4>). The quantities aT,e m = J dD,Y^ m (n)T(n) , are coeffi- 
cients of the SH expansion, and dSl = d0d<j)sm6. Under reflection (or P) symmetry 
(n —> — n), these coefficients behave as ax,im ~> (~ 1) o>T,em- CMB physics does 
not distinguish between even and odd multipoles^3E5Zl. Therefore the power con- 
tained in even and odd multipoles must be statistically the same. We thus define 
the quantity: 

C x 1 V C x (8) 

where C x are power spectral estimates for X = TT, TE, EE and BB. The sum is 
meant only over the even or odd I and this is represented respectively by the symbol 
+ or — . Therefore, two estimators can be built from Eq. the "ratio" R x = 
C x /C x (see [1361 [I37l EH]) and the "difference" D x = C X -C X (see [1381 EQ]), 
where C± is the band power average contained in the even (+) or odd (-) multipoles. 
In Fig. [7] (reprinted from |135| ) we plot the percentage related to the WMAP 7 year 
P anomaly for TT versus £ m ax in the range 10—40 for the two considered estimators. 
As evident there is not a single £ max for which the TT anomaly shows up, but rather 
a characteristic scale in the t range [15, 25]. We confirm the previously reported P 
anomaly in TT in the range A£ = [2, 22] at > 99.5% C.L.. Planck will not improve 
the signal-to-noise ratio in this range for the TT APS, since it is already cosmic 
variance dominated in the WMAP data. However, Planck has a wider frequency 
coverage and this will improve the component separation layer in the data analysis 
pipeline. Moreover, Planck is observing the sky with a totally different scanning 
strategy and this represents a benefit for the analysis of systematic effects. 
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7.1.2. Cosmological birefringence 

Linear polarization maps are components of a rank two tensor^^ and are de- 
composed by the spin harmonics a±2,£m = 

/ dfiy ±2^m(«) ± #)), where 

Y±2,foi(^) are SH of spin 2 and a±2,e m arc the corresponding coefficients. It is 
then useful to introduce new coefficients as linear combinations of the previous: 
«E,f m = -{a2,i m +a-2,im) /2 £mda B ,e m = —(a2jm—a>-2,im)/2i- These have opposite 
behaviors under a P transformation: dE,im — > ( — 1)^ aE,imi a B,im — > (— 1)^ +1 a,B,tm- 
If P is conserved, by combining the previous transformation one immediately derives 
that the cross-correlations Cj B = (a^ i m dB,i'm') an d Cf B — (a* E em aB,e'm' ) must 
vanish. Further details can be found in [1411 1142j and explicit algebra is set forth in 
the Appendix of |138j . Parity violation could, however, may change this scenario. A 
popular model for which parity is broken in the photon sector is the Chern-Simons 
perturbation to the Maxwell LagrangiarP^ 1 : AC = - \p^ l,pa F pa A v , where F * v is 
the Maxwell tensor and A p the four-potential. One of the consequences is in vacuo 
dispersion of photons, in particular those from the CMB and the rotation of their 
polarization planes, observable through TB and EB correlations, that acquire a 
signal modulated by a (or "rotated" JHUHEHl. 

The WMAP tearrPH reported a WMAF7yr = -0.9° ± 1.4° at 68% C.L.. Our 
constraint, obtained at low resolutiort^SI and considering the same estimator that 
has been used in [HE], reads a = -1.6° ± 1.7° (3.4°) at 68% (95%) C.L. for A£ = 
2 - 47. Considering Al = 2 - 23 we obtain a = -3.0°l|jo at 68% C.L. and a = 
-3.0°tt;?° at 95% C.L.. This is the same multipole range considered by the WMAP 
team at low resolution in |147j (the only other result available in the literature at 
these large angular scales) where with a pixel based likelihood analysis they obtain 
a VfUAP7yr = _ 3 go ±5 2 ° at 68 % C.L.. In [H$] it is claimed that the improvement 
expected for the Planck satellite^ in terms of sensitivity^^ is around 15. Almost 
the same number is obtained in Gruppuso et al. (135] . Both forecasts are provided 
considering just the nominal sensitivity whereas the uncertainties coming from the 
systematic effects are not taken into account. 
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Fig. 8. Polarization map (red vectors) of 0.93° X 0.93° integrated in redshift of one sky patch 
realization for pqiCMB (left panel) and p/3|SZ (right panel). For comparison, we put the logarithm 
of Compton ^/-parameter (color scaled). 

7.2. CMB induced polarization from single scattering by clusters 
of galaxies and filaments 

We discuss here two types of secondary polarization effects arising from single scat- 
tering of the CMB photons by ionized gas. These are the CMB quadrupole induced 
polarization (pqiCMB), which couples the gas density with the CMB quadrupole 
component, and the polarization induced by the gas motion transverse to the line 
of sight (pft 2 SZ). 

Instead of using an analytical model adopted in many previous studies or sim- 
ulations of individual galaxy clusters (see e.g. [1511 11521 11531 1154) ). we used high- 
resolution N-body/Hydrodynamic simulations featuring adiabatic gas physics and 
a novel box-stacking scheme that allows to reconstruct the CMB quadruple com- 
ponent and the physical properties of the scattering media along the light cone 
traversed by radiation. We generated 28 random sky patches integrated along the 
light-cone, each of about 0.86 deg 2 at angular resolution of 6" . The primordial CMB 
quadrupole information in each simulation box is computed by inverse Fourier trans- 
form of CMB quadrupole components in Fourier space at all conformal times re- 
quired by the map-making strategy. For each individual Fourier mode, we decide the 
initial value (k) by drawing a random number from a Gaussian distribution with 
variance given by the initial power spectrum obeys a power law, P*(fc) = Ak ns ~ A 
with A a normalization factor and n s a spectral index of the scalar perturbations. 
The time evolution of CMB quadrupole in each individual mode is computed using 
the CMBFastP 11 Boltzmann code (see [156] for more details). Here we focus on the 
characterization of the polarized signals in the simulations and the study of their 
statistical properties at high angular resolution. 

The results from the pixel distribution in the frequency independent maps, show 
that the linear polarization degrees follows, in logarithm scale, nearly Gaussian dis- 
tributions, centered around 10 -8 and 10 -10 for pqiCMB and p/3 t 2 SZ, respectively. 
Our simulations confirm that the polarization degree of the pqiCMB is a close 
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proxy of the electron density column and that the polarization angles of this effect 
are closely aligned due to the slow variation of primary CMB quadrupole on small 
sky patches and along the line of sight (see Fig. [sj reprinted from |156j , left panel) . 
The effect of the gas overlap along density columns causes galaxy clusters and other 
bound objects to be less prominent (with respect to the mean background) than in 
the case of the thermal SZ effect, where the signal in clusters is boosted by the high 
temperature. In the case of p/3 t 2 SZ (right panel of Fig. |8|, the polarization degrees 
and angles are weighted by the transverse velocity of the scattering media, there- 
fore the integration along the line of sight can erase contributions from individual 
collapsed objects, depending on their internal velocity structure and the effect of 
gas overlap. 

By producing maps of these secondary induced polarization effects at different 
frequencies, we confirm the strong dependence on frequency of both signals, espe- 
cially in the case of p/3 2 SZ, for which the mean value increases by a factor of ~ 100 
from the 30 GHz to 675 GHz. The high magnitudes of both signals at high frequen- 
cies may allow its detection with the next generation of sub-millimeter instruments. 

The redshift distribution of the polarization degrees shows that the contribution 
for the polarization signal is highest at z ~ 1 and z ~ 0.5 for the pqiCMB and 
p/3 2 SZ, respectively. Finally, only about 7% of the total signal comes from z > 4 for 
the former and z > 3 for the latter and both signals converge rapidly at larger z. 

7.3. Mixing of blackbodies: creation of entropy and dissipation of 
sound waves in the early Universe 

There is a very important connection between the spectrum of the monopole or 
sky averaged CMB, which is an almost perfect blackbody and COBE/FIRAS 1 ^ 
detected no deviation from Planck spectrum, and the angular anisotropics precisely 
measured by WMAP^, SPT^l, ACT^SHl, Planck and other experiments on scales 
corresponding to comoving wavenumber 10~ 4 <J k <J 0.2 Mpc~ , including the 
damping tail due to photon diff usioriMULMHI . The power that disappears from the 
CMB APS because of Silk damping appears in the energy spectrum of monopole as 
jp^H, yU and intermediate- type distortion£HEHl. The primordial power spectrum, 
at comoving wavenumbers 8 <J k <J 10 4 Mpc -1 (mostly inaccessible by any other 
means), can thus be recovered by precise measurements of the energy spectrum of 
the monopole. 

The current constraints on the primordial power spectrum, including Ly-a forest 
constraintSGM are shown in the left panel of Fig. [9] (reprinted from |165j ). At 
present the small-scale constraints from COBE/FIRAS y-type (2<r limit y <J 1.5 x 
10 -5 ) and n-type (u <J 9 X 10~ 4 ) distortions are very weak and considerable freedom 
is allowed on small scales. Proposed future experiment PIXIEP^ would improve the 
small-scale constraints by a factor of ~ 2500 and start probing the interesting region 
of the parameter space, extending our knowledge of the primordial power spectrum 
by many orders of magnitude in terms of the scales probed. 
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Fig. 9. Left panel: current and future constraints on primordial power spectrum as a function 
of comoving wavenumber k. Possible spectra on small scales allowed by current data are also 
shown. Right panel: the spectrum resulting from mixing of blackbodies. We have used the linear 
(in (AT/T) 2 ) solution to make the plots but used a large value of AT/T to make the differences 
visible. Effective temperature of the spectrum defined by writing the occupation number as n = 
l/(e hv i (^B^off) — 1) as a function of dimensionless frequency, x = hu/kT is plotted. At high 
redshifts, z > 10 s , the spectrum comptonizes rapidly to create a fi-type distortion or Bose-Einstcin 
spectrum, also shown in the figure. 



Previous calculations of spectral distortions in CMB from Silk damping^2H121J 
underestimated the energy in sound waves and also assumed that all the dissipated 
energy goes into creating spectral distortions. The physics of creation of spectral dis- 
tortion becomes very simple if we consider the fact that diffusion of photons, which 
damps the CMB perturbations, is in fact mixing blackbodies of different temper- 
atunGHEZH. The right panel in Fig. [9] (reprinted from [173] ) shows the result of 
averaging two blackbodies with temperatures T ± AT. The resulting spectrum is 
marked 'Average(Y)' and is a y-type distortiorP^ on top of a blackbody with tem- 



perature T 



1 + {AT/ry 



with the two curves crossing at x = 3.83. The averaging 
of two blackbodies adds energy as well as photons to the average CMB monopole 
and therefore not all the energy can be used to create spectral distortions. It is 
straightforward to show, by using Taylor series expansion of the initial blackbodies 
up to second order in AT/T and then doing the ensemble average, that only 1/3 
of the dissipated energy goes into y//Lt-type distortions^^ and 2/3 just raises the 
temperature. Applying the above procedure to CMB immediately gives us the rate 
of energy injection into CMB and the resulting /i distortiort^l, 



d/i 

dt 



1.4 



distortion 



d AE 



distortion 



-4k 



d 

—2.8 

dt 



k 2 dk 
k 2 dk 



P l (k) 



1=0 



f^PiW [©0 + 36?] 



(9) 



where is the fractional energy going into CMB distortion, 0£ are the multipole 
moments of CMB temperature perturbation transfer function, P t is the initial power 
spectrum, and in the second line we have used the fact that during tight coupling 
the £ > 2 modes are suppressed. The time derivatives are easily calculated using the 
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tight coupling solutions with Silk damping or by using the first order Boltzmann 
equation. A nice feature of the approach presented above is that the energy injected 
into the distortion can be directly identified with the increase in entropy of CMEP^Hl. 
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